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^Integral  scales  for  fully  developed  turbulent  flows  may  be  defined  as 
moments  of  the  two-point  tensor  correlations  among  the  dynamical  variables 
which  characterize  the  flow,  for  example  velocity  and  temperature.  We  show 
In  this  paper  that  this  concept  of  Integral  scales  allows  for  a  rigorous 
derivation  of  rate  equations  for  the  integral  scales  (at  a  single  point) 
if  two  major  assumptions  are  made.  The  first  assumption  Is  that  the 
Integral  scales  are  sufficiently  small  (compared  with  the  scale  of  spatial 
variations  of  the  mean  quantities)  to  enable  one  to  neglect  the  difference^. 
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^between  convection  at  two  points  separated  by  less  than  the  Integral  scale 
Itself.  This  assumption  Is  necessary  to  obtain  convections  of  the  Integral 
scale  (which  Is  a  one-point  function)  at  a  single  point  and  thus  represents 
the  most  general  physical  requirement  for  the  very  exlstenpe  of  closed 
equations  for  the  Integral  scales.,  The  assumption  which  Is  mathematically 
embodied  by  the  moment  expansion  of  the  correlation  tensor  4 Section  1) 
restricts  the  rigorous  theory  to  nearly  homogeneous  flows.  We  shall  see 
(Section  2)  that  our  theory  Is  In] satisfactory  agreement  with  the 
experimental  evidence  for  nearly  homogeneous  flows.  We  follow  convention 
in  making  the  bold  suggestion  thai  the  theory  be  applied  outside  the 
rigorous  limits  of  Its  derivation!  for  at  least  the  reason  that  at  the 
present  there  seems  to  be  no  sound  theoretical  alternative  to  solving  for 
the  full  two- point  tensors  jlnjjpia  11  when  near- homogeneity  falls. 

v — >The  second  major  assumption  that  we  find  necessary  for  a  rigorous 
derivation  of  rate  equations  for  the  Integral  scales  stems  from  a  coupling 
between  the  evolution  of  the  Reynolds  stress  and  other  one-point 
correlations^)  the  evolution  of  the  Integral  scale..  The  nature  of  this 
coupllnayfi,  In  our  opinion  /only  partially  understood  at  present.  To 
proceed J9«e  observe  that  we  can  obtain  from  the  two-point  correlations  at 
every  spate  point  and  at  every  Instant  of  time,  both  the  Reynolds  stress 
(by  collapsing  the  two  points  to  one),  and  the  Integral  scale  (by  a 
suitable  Integration  In  relative  separation).. 

In  view  of  this  fact.  It  appears  natural) to  propose  that  the  concepts 
of  second-order  closure,  which  have  been  successfully  developed  over  the 
past  decade,  be  extended  from  applying  to  one-point  correlations  to  apply 
generally  to  the  two-point  correlations.  It  Is  understood  that  the 
"simplest11  proven  equations  are  to  be  considered  for  the  relevant  models 
until  these  are  proven  unsatisfactory.  We  shall  see  In  what  follows  that 
this  assumption  offers  a  valuable  guide  to  obtaining  rate  equations  for  the 
Integral  scale  and  allows  us  to  Incorporate  In  the  generalized  theory  all 
the  major  successes  of  currently  available  second-order  closure. 

It  may  be  worthwhile  to  point  out  that  we  do  not  derive  the  scale 
equations  from  the  Reynolds  stress  equations.  Rather  both  equations  follow 
from  an  equation  for  the  two-point  correlation  tensor. 
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PREFACE 


There  is  more  than  one  way  to  look  at  any  physical  problem.  In  the 
field  of  turbulence  research,  it  has  become  fashionable  to  define  almost  any 
unsteady  flow  in  the  wake  of  an  unaccelerated  body  as  turbulence.  A  lot  of 
these  unsteady  flows  contain  large-scale,  coherent  structures  and,  thus, 
there  has  arisen  a  "large-scale  eddy  cult."  The  basis  thesis  of  this  group 
is  that  turbulence  should  or  (for  the  extremists)  must  be  attacked  by  some 
technique  which  identifies  the  turbulence  with  the  interaction  and  decay  of 
such  large-scale  structures.  A  corollary  of  this  position  is  that  closure 
methods  do  not  or  (for  the  extremists)  cannot  address  themselves  specifi¬ 
cally  to  large-scale  eddies  and,  therefore,  are  not  really  anything  but 
dull,  unphysical,  and  temporary  methods  for  dealing  with  turbulent  flows. 

The  authors  of  this  paper  do  not  believe  that  these  people  really 
understand  the  nature  of  closure  calculations  at  the  present  time.  Not  only 
have  closure  methods  demonstrated  the  existence  of  large-scale  eddies  in  two 
cases  (the  roll  eddies  of  the  marine  planetary  boundary  layer  and  unsteady 
large  eddies  in  the  flow  behind  a  rearward-facing  step  (see  references  11 
and  12  in  the  main  body  of  this  report)),  but  these  eddies  have  been  resolved 
in  all  their  gory  detail.  This  can  be  accomplished  when  the  closure  equat¬ 
ions  are  used  in  their  elliptic,  time-dependent  form  and  the  grid  spacing  is 
fine  enough.  Why  is  this  so?  It  is  because  the  Euler  equations  (which 
govern  the  formulation  and  a  great  deal  of  the  behavior  of  large  eddies)  are 
contained  in  the  time-dependent,  elliptic  equations.  Thus  we  submit  that 
closure  techniques  not  only  can  describe  large-scale  eddies  but  must  do  so 
if  the  time-dependent,  elliptic  forms  of  the  equations  are  used  and  the  grid 
spacing  is  small  enough  to  resolve  these  eddies. 

One  further  point.  The  Karman  vortex  wake  is  an  unsteady  flow  associa¬ 
ted  with  an  unaccelerated  body.  We  do  not  prefer  to  think  of  it  as  a  form 
of  structured  turbulence  (nor  did  those  who  first  studied  the  phenomenon) 
although,  if  one  were  a  member  of  "the  cult,"  this  point  of  view  might  be 
taken.  The  reason  that  we  choose  not  to  consider  it  as  structured  turu- 
lence  is  that  there  is  ample  evidence  of  its  existence  in  laminar  flow. 
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From  the  closure  point  of  view,  the  Kerman  wake  Is  looked  at  as  an  unsteady 
flow  peculiar  to  the  body  that  produces  It.  This  unsteady  flow  Interacts 
with  Itself  to  decay  either  through  laminar  exchange  and  dissipation  or 
through  turbulent  exchange  and  dissipation,  depending  on  the  Reynolds  number. 

The  research  reported  here  Is  a  description  of  our  first  attempts  to  make 
a  closure  theory  of  turbulence  that  Is  compatible  with  the  large-scale  structures 
that  we  must  Inevitably  find  when  we  run  our  closure  codes  In  an  elliptic 
manner.  The  work  reported  Is  one  completed  step  In  this  direction.  As  noted 
later  In  the  text,  we  do  not  consider  this  work  complete.  We  sincerely  regret 
that  we  had  to  terminate  this  work  before  our  attempt  to  construct  a  more  general 
closure  formulation  could  be  completed.  While  we  consider  many  of  the  detailed 
large-eddy  studies  to  be  outstanding,  we  also  believe  that  elliptic,  non¬ 
steady,  closure  techniques  will  be  the  backbone  of  turbulence  computations 
that  will  be  of  use  to  the  military  for  the  next  twenty  years,  and  that  closure 
techniques  that  are  compatible  with  this  approach  should  and  will  be  pursued. 
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1.  INTRODUCTION 


Thinking  of  the  Integral  scales  of  turbulence  as  moments  of  the 
two-point  correlation  tensors  generates  a  general  approach  to  the  dynamical 
determination  of  these  scales.  The  purpose  of  this  paper  is  to  demonstrate 
the  validity  of  this  proposition  and  to  develop  a  number  of  its 
consequences. 

o 

In  order  to  calculate  the  level  of  turbulent  kinetic  energy,  q  /2,  the 
separate  levels  in  the  "energy  components,"  ( u* 2 ,  v'2  and  w'2),  as  well  as 
the  turbulent  stresses  (e.g.,  u V ) ,  the  second-order-closure  approach  has 
been  developed  to  a  high  degree  over  the  last  decade.  To  make  calculations 
of  the  turbulent  fluctuations  indicated  above,  it  is  necessary  to  provide 
information  on  the  behavior  of  the  turbulence  scales  that  have  been 
introduced  in  second-order  models  that  represent  turbulent  transport, 
isotropization,  and  dissipation.  It  is  desirable  to  obtain  local  rate 
equations  for  the  turbulence  scales  because,  for  example,  one  of  them 
represents  the  size  of  typical  energy-containing  eddies  and  this  size  varies 
considerably  from  point  to  point  in  a  turbulent  flow.  A  similar  observation 
applies  to  the  size  of  the  eddies  that  are  mainly  responsible  for  the 
dissipation  of  the  turbulence.  Ue  develop  in  this  paper  a  technique  to 
determine  local  rate  equations  for  the  integral  as  well  as  the  microscales 
of  turbulent  flows.  We  follow  and  generalize  the  basic  ideas  introduced  by 
6.  I.  Taylor1  who  focussed  his  attention  on  the  important  case  of  isotropic 
turbulence.  The  normalized  autocorrelation  function  f(r)  of  the  component 
of  the  velocity  In  the  direction  £with  the  same  velocity  component  at  a 
distance  r  has  schematically  the  form  indicated  In  Figure  1.  Here  f(r)  is 
defined  to  be  the  normalized  time  (or  space,  or  ensemble)  average 

(u|.{x)  u’<: 
f(r)  =  .  — ~  r:  T— 

(4Z(x 

The  typical  curve  Is  not  a  Gaussian  except  when  the  turbulence  is  very  weak 
(e.g.,  In  the  final  stages  of  decay  of  Isotropic  turbulence).  It  Is 
therefore  necessary  to  allow  for  two  distinct  spatial  scales  n  order  to 


9 


t  \|  \  O l'L  1  3LA 


ly  characterize  turbulent  flows,  and  this  was  indeed  done  by 
le  smaller  scale  ("microscale")  is  defined  by  the  curvature 
;he  origin 


\2  =  lim  C-f"( r)]”1  >.0 
r*o 

shown  to  be  characteristic  of  the  eddies  responsible  for  viscous 
i.  The  longer  scale  ("integral  scale")  is  defined  by  Taylor  in 
as  to  account  for  the  non-Gaussian  tail  of  f  as  the  area  under 
t,  namely. 


L 


f  = 


f(r)  dr 


lies  degenerate  into  a  single  scale  as  f  becomes  Gaussian.  The 
i  between  the  two  scales  can  also  be  inferred  from  the  form  of  the 
rum  of  the  turbulent  field  because  the  power  spectrum  is 
easily  related  to  the  Fourier  transform  of  f(r).  One  geometric 
the  two  scales  \  and  Lf  is  shown  in  Figure  1  where  Jl  \  is  the 
f  the  parabola  tangent  to  f  at  the  origin  with  the  r  axis,  and  If 
ructed  that  the  hatched  areas  (denoted  by  A)  are  equal.  We  shall 
ylor's  definitions  yield,  when  appropriately  generalized, 
which  transform  as  tensors  and  are  therefore  meaningful  in  flows 
general  geometry.  We  shall  then  show  that  the  scales  satisfy 
d  dynamical  equations  which  we  shall  derive  from  the 
es  equations  for  fluid  flow. 

pose  a  generalized  second-order  closure  which  makes  the  equations 
■point  correlation  functions  self-contained.  Thus,  for  the 
ilocity  correlation  tensor  R^j,  we  say  that  we  have  a 
*  closure  if  we  can  write 
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aRij  -  r  rR  I 
Tij  LRkA-1 


where  T  is  a  tensor  functional  of  R^j.  We  shall  restrict  T  considerably  by 
requiring  integrals  to  appear  only  through  the  integral  scale  tensor 
which  we  define,  following  Rotta2,  as  a  weighted  moment  of  Ru 


"3“  Aij^c^  *  f  $rir2  ^-c* 


The  two  variables  xc  and  r,  are  defined  below.  We  observe  that,  for 
isotropic  turbulence 


(is) 

Aij 


(is) 

Akk 

3 


The  normalization  factor  has  been  adopted  so  that,  for  isotropic  turbulence 


.  _  1  .  s  Lf  *  2Lg  s  2  , 

A  =  3  Akk  "  3  ~  3  Lf 


It  is  interesting  to  note  that  this  choice  of  normalization  leads  to  a 
scalar  scale  a,  which  is  2/3  of  the  longitudinal  integral  scale  and  that  the 
scale  now  defined  is  approximately  equal  to  that  generally  used  in  current 
second-order-closure  calculations.  The  basic  reason  for  the  definition 
adopted  above  for  (generalized)  second-order  closure  can  now  be  made  clear. 
The  tensor  R^j  contains  the  information  needed  to  obtain  the  kinematic 
Reynolds  stress 


u}ui  *  lim  R^j  U,£) 

and  also  the  scales  that  enter  the  modeled  terms  in  the  equation  for  u{uj  if 
these  are  assumed  to  be  related  to  the  Ajj  just  introduced. 
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In  order  to  derive  local  rate  equations  for  A^j,  we  find  it  necessary 
to  generalize  the  standard  Integral  methods  to  extract  average  information 
from  equations  in  many  independent  variables.  The  procedure  that  we 
introduce  Is  shown  to  coincide  with  the  standard  integral  methods  where 
these  are  applicable.  Our  generalization  is  made  necessary  in  order  to 
treat  the  two-point  correlation  tensors  of  the  theory  of  turbulence  because 
these  are  convected  independently  at  two  distinct  points.  The  use  of  the 
moment  expansion  (which  is  necessary  to  obtain  rate  equations  for  scales  at 
a  single  point)  restricts  the  ensuing  theory  to  nearly  homogeneous  flows. 

We  show  in  the  following  that  our  theory  is  in  satisfactory  agreement  with 
the  experiments  on  nearly  homogeneous  flows.  We  believe  that  our  theory 
(like  the  Chapman-Enskog  theory  of  molecular  transport)  will  give  useful 
information  outside  the  limits  of  its  rigorous  derivation. 

1.1  Equations  for  the  Two-Point  Correlations 

We  consider  incompressible  flows  with  constant  density  which  may  carry 
a  passive  additive.  We  assume  that  the  governing  equations  are  the  Navier- 
Stokes  equations 


-2L. 

dx^ 


ui 


the  continuity  equation 


(1) 


&ui 

ax.j 


o 


(2) 


and  the  convection-diffusion  equation 


at 


(3) 


The  notation  adopted  is  that  of  cartesian  tensors  with  the  following 
variables:  the  eulerlan  velocity  Is  u^,  a  field  function  of  the  Independent 
variables  x  *  (x^)  and  t;  the  kinematic  pressure  is  denoted  by  p  which  is 
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the  ratio  of  pressure  to  mass  density;  the  passive  additive  is  measured  by 
9  which  represents,  for  example,  the  temperature.  The  two  transport 
coefficients  are  v,  the  kinematic  viscosity,  and  0,  the  diffusivity  of  the 
passive  additive.  The  Laplacian  operator  at  the  point  x_  is  written  as  vj* 
and  it  is  defined  by 


2 


V 


X 


(4) 


We  introduce  the  Reynolds  decomposition  of  the  dependent  variables  by 
writing 


-  I 

U.J  =  Uj  +  Uj 


(5) 


P  *  P  +  P' 


(6) 


0*9+9'  (7) 

where  the  average  can  be  understood  as  an  ensemble  average  or  as  a  space  (or 
time)  average.  In  this  latter  case  it  is  assumed  that  there  is  a  marked 
separation  between  the  short  (or  small)  scale  of  the  turbulent  fluctuations 
and  the  large  (long)  scale  of  the  variation  of  mean  quantities.  The 
averages  are  thereby  always  interchangeable  with  the  space  and  time 
derivatives.  The  transport  coefficients  are  assumed  constant.  For  any 
arbitrary  quantity  A,  we  shall  use  Interchangeably  the  overbar  or  the 
bracket  to  denote  the  turbulent  average.  Thus 

A  =  <A>  (3) 


We  now  Introduce  the  following  two- point  correlation  tensors 

R(x,  i)  *  <0'(x)  0'(£)>  =  R  (9) 


14 


R^(x,  £)  *  <u}(x)  e'(^)>  s  R^ 


(10) 


Rij  i)  *  ^(x)  =  Rij 


(11) 


The  order  of  Indices  and  arguments  is  essential  because  the  correlation 
tensors  are  not  symmetric  under  separate  exchange  of  indices  and  of 
positions.  Note  that  when  differentiating  the  argument  can  be  omitted 
without  ambiguity  when  a  function  of  a  single  argument  is  differentiated. 
Note  also  that  the  Reynolds  stress  <u}uj>  is  the  limit  of  R^j  as  £  ♦  x. 

Substituting  the  Reynolds  decompositions.  Eqs.  (S)  -  (7),  into  the 
basic  Eqs.  (1)  -  (3)  we  obtain,  for  the  momentum  equation 


dui  aui  _  aut  -  aui 
at  at  k  axk  k  axk 


+ 


-2L  +  3£l 

dXj  axi 


ui  + 


V  ul 


for  the  continuity  equation: 


au<  au- 
— I  +  — L  x  o 
axj  ax^ 


and  for  the  convection-diffusion  equation: 


(12) 


(13) 


M  +  ael  +  -  _ae.  + 
at  at  k  axk 


uk 


89 1 

axk 


DV*  0  ♦  D  vj  e' 


(14) 


We  average  the  "total"  Eqs.  (12)  through  (14)  and  remember  that  the 
averages  are  so  defined  as  to  commute  with  the  space  and  time  derivative. 
The  result  in  "mean  convection"  form  is,  for  momentum 
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M. 

dxi 
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<uuui> 


+ 
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For  the  continuity  condition  we  have 


(15) 


au1 

ax7 


o 


(16) 


whlcn  was  used  to  give  the  form  exhibited  for  the  right-hand  side  of  Eq. 
(15).  For  the  passive  additive  we  find 


ae 

at 


+ 


uk 
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axk 


—  <ui  e*>  +  D  v?  e 
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We  remember  Reynold's  basic  observation  that  turbulence  gives  effective 
transport  effects  to  the  fluid  in  the  form  of  a  kinematic-stress  tensor, 
<Uy  u|>,  and  of  a  flux  vector,  <u^  e‘>,  for  the  passive  additive. 

Subtracting  the  mean  from  the  total  equations,  we  obtain  the  starting 
equations  for  the  determination  of  the  rate  equations  of  the  correlation 
tensors  given  by  Eqs.  (9)  through  (11).  These  equations  are  the 
fluctuation  equations: 


for  momentum 


auj  _ 
— -  +  uk 
at  K 


axk  dx^ 


auj 


-  “i  ^  ^ +  v  v' 


X  ui 


(18) 


for  continuity 


ax1 
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and  for  the  fluctuations  of  the  passive  additive 


(20) 


58*  a 

axk  ’  axk 


[u^e1  - 


<ujce'>]  +  o  vj-e' 


The  fluctuation  equations  have  been  put  in  the  "mean  convection"  form  in 
order  to  emphasize  that,  for  the  fluctuations,  the  governing  equations  are 
of  the  primitive  form  with  two  modifications:  turbulent  transport  tensors, 
[Uy  u;  -  <Uy  u)>]  and  [u^  e‘  -  <u^  e‘>],  respectively  and  also  turbulent 
"production"  -u^  (auij/dxu)  and  -u^  {ai/axu),  respectively. 

We  can  now  derive  rate  and  continuity  equations  for  R^j.  Multiply  the 
au}(x)/at,  Eq.  (18),  by  Uj(^)  and  the  equation  for  aUj(^)/at  (derived  from 
Eq.  18)  by  u{(xh  Then  add  and  average.  Using  the  definition,  Eq.  (11), 
we  obtain 


(21) 


We  observe  that  the  index  1  systematically  accompanies  the  argument  £  and 
that  the  index  j  systematically  accompanies  the  argument  The  terms  that 
require  modeling  are  those  for  which  does  not  appear  explicitly.  They 
are  the  pressure  velocity  correlation  and  the  triple  velocity  correlation. 
It  is  of  interest  to  note  a  first  important  difference  between  the  familiar 
equation  for  the  Reynolds  stress  tensor  <u}uj>  and  Eq.  (21)  for  the 
two-point  correlation  tensor  R^j.  In  the  equation  for  the  Reynolds  stress, 
the  viscous  term  requires  closure  since  this  is  given  by 


Aim  L  (yj  +  72)  R  U,£) 

£*x  L 


(22) 


which  cannot  be  expressed,  in  general,  in  terms  of  u}uj  alone.  By  contrast, 
the  viscous  term  in  the  equation  for  R^  does  not  require  modeling. 
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In  addition  to  the  momentum  equation,  Eq.  (21),  R^j  also  satisfies  two 
continuity  equations.  Thus  multiplying  Eq.  (19)  by  Uj(^)  we  have 


6ui(x) 

—  =-  =  0 


ax 


i 


which,  when  averaged,  reads 


(23) 


(24) 


or 


RU  <*’  *>  *  0 

Similarly 

— R44  (x,  y )  =  0  (26) 

ayj  n  - 

Equations  (25)  and  (26)  constitute  a  second  and,  in  our  view,  very  Important 
difference  between  the  equation  for  the  one-point  Reynolds  stress  tensor  and 
the  two-point  R^j.  In  fact  the  stress  ujui  does  not  satisfy  a  continuity 
equation  as  Is  well  known  from  Llghthill's^  approach  to  the  problem  of 
determining  the  noise  generated  by  turbulence.  By  contrast,  R^j  satisfies 
Eqs.  (25)  and  (26)  which  will  be  seen  to  be  stringent  requirements  on  the 
modeling  of  the  R^j  rate  equation. 

By  a  reasoning  entirely  parallel  to  the  one  outlined  for  the  derivation 
of  the  rate  equation  for  R^j,  we  can  obtain  the  rate  equation  for  the  flux 
vector  Rj  of  the  passive  additive.  The  result  of  the  calculation  is 

*  [“*<*>  £  *  “*'*>  £]  Ri‘i*  0  *■<*.<>  ■ 
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(27) 


1/ 

I 


du< 


Rik  r1 '  Rk(i*  p 

1K  dyk  K  axk 


-  ~  <uk(x)uj(x)  e'{£)> 


—  <u}(x)uj^>  e' {%)> 


+  (v  vj  +  D  V*)  Ri(x,  ^) 


The  reader  can  convince  himself  that  no  simple  notational  convention 
unambigously  suppresses  the  arguments  of  the  flux  for  the  production  term; 
hence  the  apparent  pedantry  in  carrying  the  arguments  of  explicitly. 

This  result  also  applies  to  the  continuity  equation  for  the  flux  which  reads 


Ri(x,  y)  =  0 
axi  1  *- 


(28) 


We  observed,  when  discussing  Rjj,  that  the  viscous  term  and  the  continuity 
conditions  make  for  substantial  differences  between  the  equation  for  u}uj 
and  for  R^j.  The  same  concepts  apply  to  Eqs.  (27)  and  (28):  (1)  The 

viscous  term  in  the  equation  for  Rj  is  automatically  closed  and,  (11)  the 
two-point  flux  equation  must  be  compatible  with  the  continuity  requirements, 
Eq.  (28).  The  physical  meaning  of  the  continuity  requirements  on  R^j  and  on 
Rj  Is  best  seen  by  differentiating  Eq.  (27)  with  respect  to  d/dXj  and  using 
Eq.  (28).  We  obtain 


0  + 


^  Mi. 


_ 2i_ 

ax^axfc 


cuj^xJujOOe'  (^)> 


+  0  - 


aRkU»j£)  5Uj 
axj  axj^ 


(29) 


Equation  (29)  rearranges  to 


5u<  5Rk(x,y)  .2 
H  <p'(*)e'(^)>  -  -2  —  — ~ - (30) 


8Xk  dX1 


ax1axk 


which  can  be  compared  with  the  familiar  Poisson  equation  for  the  pressure 


19 


fl  uctuatlons 


5u^  auj  dUf  dUj  d2<UfUj> 

1  axj  ax^  axj  ax^  ax^axj 


obtained  from  the  momentum  equation  for  fluctuations,  Eq.  (18),  by  using  the 
Incompressibility  of  velocity  fluctuation,  Eq.  (19).  In  fact,  Eq.  (31) 
multiplied  by  e‘ {^)  and  averaged  yields  Eq.  (30).  An  analogous  reasoning 
applies  to  Rjj.  In  this  case  using  Eq.  (21)  together  with  the  continuity 
condition  Eq.  (25)  gives 


v2  <p'(x)u^{£)> 


2  aux  aRpk 


a2 

ax^axc 


<U^(_x)Up(x)u^(j^)> 


This  equation  can  also  be  obtained  by  multiplying  Eq.  (31)  by  uj^(j^)  and 
averaging.  Vie  thus  see  that  Eqs.  (30)  and  (32)  contain  essential 
Information  on  the  long-range  properties  of  the  pressure  fluctuations  as 
given  In  detail  by  the  Poisson  equation  (31).  We  believe  that  a  model  for 
Rjj  that  successfully  Incorporates  Eqs.  (30)  and  (31)  would  help  clarify  the 
structural  as  well  as  spectral  properties  of  turbulent  flows. 

We  complete  this  subsection  by  obtaining  the  rate  equation  for  the 
temperature  autocorrelation  function  R  =  <0*  (x)0‘ (£)>  as  defined  in  Eq. 

(10).  We  multiply  the  equation  for  ae'(x)/dt  by  e'(£)»  add  to  it  the 
equation  for  ae'(^)/at  multiplied  by  e'U)  and  add.  The  result  Is  the 
equation 

ft  *  [“k(i'  £  *  vn  R  ■  -  [Mi*'  £ +  R*(^> 

-  -S-  <ui ( x)  e'(x)0'(y)> 

axk  k  -  “ 

-  <u&(£)  e‘(x)  e’(^)> 

♦  o  (v£  +  v|)R  (33) 


We  ask  patience  for  calling  attention  again  to  the  arguments  of  the 
two-point  flux  vector.  In  parallel  with  and  Rj,  and  by  contrast  with 
the  equation  for  the  one-point  temperature  variance,  there  is  no  need  to 
model  the  viscous  term  in  the  two-point  temperature  correlation  equation. 

No  continuity  requirement  arises  for  R  since  the  velocity  does  not  enter  its 
definition. 

In  summary,  the  basic  equations  for  the  correlation  tensors  and  their 
continuity  conditions  are  given  by  Eqs.  (21),  (25),  and  (26)  for  R^j,  Eqs. 
(27)  and  (28)  for  R^  and  Eq.  (33)  for  R. 

We  consider  now  the  limit  process  that  makes  the  two  points  of  the 
correlation  function  collapse  to  a  single  point. 

For  this  purpose  we  introduce  a  convenient  coordinate  system.  The 
centroid  vector  xc  and  the  relative  position  vector  £  are  defined  by 

xc  *  j  (x  +  ^)  (34) 


£  =  £  -  x  (35) 

These  equations  are  inverted  by 

*.a25c"jl*£*-c  +  7l-  (36) 

The  geometry  of  the  transformation  is  shown  in  Figure  2.  With  these 
coordinates  the  collapse  of  the  two  distinct  points  x  and  £  can  be 
accomplished  by  taking  the  limit  £  -*■  0  holding  xc  fixed;  thus  the  two 
points  £  and  £  collapse  to  their  centroid.  To  express  the  derivatives  that 
appear  in  Eqs.  (21),  (27),  and  (33)  for  the  correlation  tensors  and  the 
continuity  equations,  (25),  (26)  and  (28),  we  make  use  of  the  chain  rules 


j_s  l  a  +  _a_ 

ay1  2  &xc1  ari 

We  note  that  for  functions,  FU)  of  x,  only,  l.e.,  functions  that  are 
Independent  of  £,  we  have 


3F(x) 

ax,- 


(39) 


Similarly,  for  functions  G(^),  which  are  independent  of  we  have 


dG(j^) 

r 


(40) 


To  illustrate  the  collapse  process,  we  consider  the  divergence  of  the  triple 
velocity  correlation  which  has  two  terms  in  the  equation  and  but  a 
single  term  in  the  equation  for  <u|uj>.  We  note  first  that  the  relation 
that  holds  in  the  collapse  limit  as  defined  above  for  the  Reynolds  tensor  is 

<u;(xc)uj(xc)>  =  Aim  R^Xj.  -  \  r,  ^  j  rj  (41) 

Further,  the  turbulent  kinetic  energy  is  defined  as  q2/2  where 

q2(xc)  =  (42) 

The  two  terms  in  Eq.  (21)  that  contain  triple  velocity  correlations  are: 


T1jt**> 


^  <ui(x)uj(^)u^(x)> 


<u}(x)  Uj(^)  u£(*)> 


(43) 


Introducing  centroid  and  relative  variables  and  using  continuity  of  velocity 
fluctuations  we  have 


7  I*  *c  + 
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=  “4  -  ?  r)]4c '  *  L)  +  ? r)/ 

*  - ? 1)  “i(ic  *  1 1)[^  “j(ic  *  ?  t)]^> 

We  can  now  readily  compute  the  limit  as  £  -*■  0  .  Thus 

Tij(*c  -  ?£•  £e  *  ?l) 


‘  Tij  <V  ic> 


\[^7ui<-c'] “*  (ic’  “j(i 
=  ^  uj(2Sc>  u*(&)) 


'j'icV  (“i<*c>  ui(Sc'[^“j«5c»]^) 


Applying  the  limit  r  +  0  to  Eq.  (21),  we  find  the  well-known  equation  for 
the  Reynolds  stress  tensor 


9  ll'  A  i.  9 

it  "<“j  *  “t  i^ 


“i“j  *  <^7  “j)  *  (“i  ^“) 


-r— r-  8U4  dU1 

*  -  U.U,  - ±-  +  - -  UkU^ 

1  u  axcu  Qxck  J 


‘  9^  “W  +  v<u*  7c  u?  +  <uj  7c  ui>  {46 

where  all  correlations  are  evaluated  at  the  centroid  vector.  Equation  (45) 
is  equivalent  to  the  second  term  In  the  right-hand  side  of  Eq.  (46). 

Application  of  the  collapse  process  to  the  continuity  equation,  Eq. 
(25),  gives 


ns  the  fact  that  the  Reynolds  stress  does  not  fulfill  a  simple 
ondition. 


ce  that  the  two  identities  commonly  exploited  to  introduce  the 
"tendency  towards  isotropy"  and  of  "turbulent  dissipation"  have 
jes  for  the  relevant  two-point  correlations. 

pressure  gradient  velocity  correlation  in  Eq.  45,  we  have  the 
identity 


»  obtain 
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'  <p'uj 
i  J 


>  + 


a 

dxcj 


<p'U;> 


("pressure  diffusion") 


axcj  dxci )/ 


("tendency  towards  isotropy") 


(48) 


(49) 


t  term  is  traceless  as  a  consequence  of  continuity.  For  the 
I  pressure  gradient- velocity  correlations  in  Eq.  (21)  we  can 
alogous  decomposition,  namely. 


ax 


cj 


(*) 


<u)(x)p' (^)> 

au}(x\l(" 
axcj  7 J 


("two- point  pressure 
diffusion") 
two-point  tendency  toward 
isotropy") 


(50) 


ition  Is  justified  by  the  fact  that  each  of  the  two  terms  in 
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square  brackets  in  Eq.  (50)  reduces  to  the  corresponding  tern  in  Eq.  (49). 
We  have,  in  fact,  for  the  two-point  tendency  towards  isotropy 


/  dui(y)  3ui(x)\ 

^'l*c  •  r-1 1^7  4; 


In  the  limit  as  r  +  0  this  becomes 


f  au'.ixj  au^x.) 

p'lx.1  ■*“£-♦ 

\  xci  ^  5xcj 


Furthermore  we  note  that  the  two-point  tendency  T^j  U,y)  is  trace  free 

Tkk><i,  =  0 


The  viscous  term  in  the  Reynolds  stress  equation  is  often  rearranged 
using  as  an  identity  Leibintz's  rule 


u}uj  =  <u}(vcUj)>  +  <Wc^)Uj>  +  2 


dUj  duk 

£xck  aXck 


which  can  be  used  to  rewrite  the  viscous  term  in  Eq.  (46)  as 


v  <U)  Vguj>  +  <Uj  7^  u-> 

=  v7^  <u)u]>  (viscous  "diffusion"  of  the  Reynolds  stress) 


\xck  3xcv 


("turbulent  dissipation") 
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In  good  analogy,  we  can  write,  for  the  two-point  viscous  term  in  Eq.  (21): 


v(vj|  +  Vy)  R^j  =  v<U-(x)  72  Uj(^)>  +  v  <(7^Uj(^))Uj(^)> 

(viscous  “diffusion"  of  R^j) 

^u)  ( x)  duj  (^) 


v  7c  Rij 


-  2v<— -  — ‘L-=-  >( two- point  turbulent  “dissipation") 

\dxck  axck  /  (55) 


The  identification  is  justified  here  as  in  the  case  of  pressure-gradient- 
velocity  correlations  on  the  ground  that  we  have  term  by  term  collapse  of 
Eq.  (55)  to  Eq.  (54). 


Before  introducing  the  models  adopted  for  our  generalized  second-order 
closure,  we  observe  that  the  collapse  of  the  two  points  in  Eq.  27  for  R^ 
gives  the  familiar  rate  equation  for  the  one-point  flux  <uje*> 
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axck 
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Similarly  Eq.  (33)  for  R  becomes  in  the  limit  r 


7c9> 


>  0: 
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dXck 
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(57) 


We  note  that  when  the  Prandtl  number  (v/D)  differs  from  unity,  the  molecular 
transport  terms  in  the  flux  equation  play  a  different  role  than  those  in  the 
Reynolds  stress  and  temperature  variance  equations. 


1.2  Model  Equations  for  the  Correlation  Tensors 
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In  order  to  achieve  our  goal  of  obtaining  dynamical  equations  for  the 
Integral  scales  associated  with  the  correlation  tensors,  we  must  close  these 
equations  In  a  suitable  form.  The  path  that  appears  most  natural  Is  that  of 
requiring  that  the  closed  equations  contain  only  the  same  two-point 
correlation  tensors  whose  Integral  scales  we  seek  to  determine.  This  point 
of  view  Is  made  plausible  by  the  successes  of  second-order-closure  modeling 
In  that  the  modeled  equations  contain  only  the  one-point  correlations  In 
question  and  scales  required  to  give  correct  dimensionality  to  the  terms 
that  require  modeling.  If  these  scales  are  interpreted  as  moments  of  the 
two-point  correlations,  no  quantity  extraneous  to  the  two-point  correlations 
needs  to  be  Introduced. 

To  specify  the  two-point  models  more  precisely,  we  adopt  two  sets  of 
rules.  The  first  set  of  four  Is  a  formalization  of  the  notion  of  invariant 
modeling  and  thus  invokes  symmetry  laws  whose  validity  is  hardly  in 
question;  the  second  set  summarizes  a  number  of  sound  dynamical  rules  whose 
validity  Is  known  to  have  limitations.  The  second  set  of  rules  is  taken 
merely  as  a  flexible  guideline,  albeit  a  very  useful  one.  The  first  set  of 
rules,  the  klnematical  rules,  is  as  follows: 

(1)  The  model  equations  are  invariant  under  space  rotations,  space 
translations  and  time  translations.  These  properties  are  of  course 
true  of  the  unmodeled  Eqs.  (21),  (27),  and  (33)  for  R^j,  Rj  and  R.  We 
thus  require  that  the  model  for  a  given  expression  have  the  same  tensor 
rank  and  symmetries  both  in  tensor  indices  and  position  labels  as  the 
expression  modeled. 

(2)  The  model  equations  are  Invariant  under  Galilean  transformations  which 
is,  of  course,  true  of  the  unmodeled  equations. 

(3)  The  modeled  equations  have  the  same  scaling  laws  as  the  unmodeled 
equations;  that  is,  the  model  for  an  expression  is  to  have  the  same 
dimensions  as  the  modeled  expression. 

(4)  The  modeled  expression  satisfies  the  same  continuity  properties  (if 
any)  that  are  satisfied  by  the  term  modeled. 
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The  second  set  of  rules  that  we  adopt  intends  to  capture  some  of  the 
broad  dynamical  features  of  turbulent  flows  extracted  by  many  workers  with 
much  labor  over  the  last  fifty  years.  They  are: 

(l1)  The  triple  correlations  represent  the  diffusive  character  of  turbulent 
flows,  an  idea  strongly  advocated  by  G.  I.  Taylor4  since  his  early  work 
on  turbulence.  There  are  likely  exceptions  to  this  rule. 

(2')  The  dissipation  occurs  at  the  smaller  scale  end  of  the  spectrum  which 
is  predominantly  isotropic.  This  idea  embodies  the  Kolmogoroff 
philosophy  and,  like  (1'),  this  rule  has  likely  exceptions. 

(3'1  The  main  effect  of  pressure  fluctuations  for  free  flows  is  to 

isotropize  the  turbulence,  the  Rotta- Satchel  or’’**  concept.  We  do 
expect  strong  exceptions  to  this  rule,  for  example  in  the  presence  of  a 
wall  or  of  strong  stratification. 

The  rules  (1'),  ( 2 ' )  and  (3')  have  proven  very  fruitful  in  second-order 
closure  of  one- point  correlations.  We  assume  them  applicable  to  the 
two-point  models  that  we  study  here. 

We  have  found  it  very  useful  to  adopt  an  additional  rule,  namely  that 
the  two-point  models  should  collapse  to  a  standard  one-point  model  when  jr  ■* 
0.  We  shall  adhere  to  this  correspondence  principle  and  adopt  as  standard 
the  second-order  modeling  of  Donaldson  and  his  coworkers  (see,  for  example. 
Ref.  7). 


We  consider  four  main  models: 
(a)  Velocity  diffusion.  We  choose 


^_<U;(x)  u£(x)  u](*)>  +-£-<u}(x)  u£(*)  u}(£): 
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*>xck 
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5RjU  +  dRJti 
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(58) 


While  this  model  slightly  generalizes  previous  models  by  allowing  for 
anisotropic  diffusion,  we  shall  be  mainly  interested  in  the  simple  special 
case 
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159) 


°kt  '  J  I'ic*  A  'ic)slce 

OjA  *  0  (60) 

where  A  Is  to  be  a  moment  of  Rjj  and  q2  *  <u|u|>.  The  quantity  r) 

Is  required  to  have  the  property 

L1j  a  0  (61) 

That  is  to  say,  L^j  does  not  contribute  to  the  evolution  of  the  Reynolds 
stress  tensor.  Thought  of  in  terms  of  the  theoretical  work  of  the  1940’s 
and  1950's,  the  term  L^j  is  associated  with  the  cascading  of  large  eddies 
into  small  ones.  In  terms  of  more  recent  work  where  both  "merging"  of 
eddies  (e.g.,  in  the  Brown-Roschko  experiments)  and  cascading  can  occur,  it 
seems  preferable  to  think  of  this  quantity  in  terms  of  eddy-size 
rearrangement.  We  shall  adopt  a  simple  linear  behavior  for  In  terms  of 
Rj j .  Thus 

•  »  J  »IJ  -  <62> 

where  V  Is  a  size  rearrangement  parameter.  Sjj  Is  required  to  satisfy 

s1j  *  u|uj  (63) 

to  guarantee  that  L^j  does  not  contribute  to  the  evolution  of  the  Reynolds 
stress  and  to  have  vanishing  leading  moment  In  order  that  L^j  give  a  growing 
contribution  to  the  Integral  scale.  This  Is  required  by  the  experiments  on 
grid  turbulence  as  we  shall  see  In  Section  2. 

(b)  Pressure  diffusion.  We  choose: 


i1 


— <p'(x)u‘i(y)>  +  — <ui(x)p‘  (y)>  = 
axci  J  axcj 


(64) 


Vie  shall  be  mainly  Interested  In  the  special  case  D3  *  0  which  seems 
adequate  at  present  mainly  because  pressure  and  velocity  diffusion  effects 
are  very  difficult  to  separate 
(c)  Tendency- towards- Isotropy.  We  choose: 


p'(x) 


auj(x) 

axci 
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(65) 


The  first  two  terms  in  the  square  bracket  are  an  obvious  generalization  of 
the  standard  Rotta  model  for  the  Reynolds  stress  rate.  The  tendency 
correction  T^j  is  required  to  have  vanishing  leading  moment  so  that  the 
entire  term  gives  a  tendency  to  isotropy  for  the  scale  tensor  and 
furthermore  T^j  must  insure  that  continuity  is  fulfilled  for  the  R^j  model. 
It  may  be  shown  that  for  isotropic  turbulence  Tjj  is  uniquely  determined, 
and  it  is  such  as  to  make  the  pressure-velocity  correlation  vector  vanish  as 
appropriate  to  this  case.  A  complete  determination  of  T^j  for  an  arbitrary 
homogeneous  turbulence  is  not  available  at  present. 

(d)  Turbulent  dissipation.  Utilizing  our  correspondence  principle  and  the 
guiding  rule  ( 2 ' ) ,  we  choose  (for  high  Reynolds  number) 


au1 ( *)  Mod 

axck  axck 


(66) 


where  b  *  1/8  (see  Ref.  7).  This  model  will  be  seen  in  Section  2  to  be 
adequate  for  near-homogeneous  shear  turbulence  at  high  (turbulent)  Reynolds 
number. 
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1.3  Moment  Expansion 

Me  propose  below  a  "moment"  expansion  for  functions  which  are  sharply 
peaked  about  a  definite  point  which  for  mathematical  discussion  can  be 
chosen  to  be  the  origin.  The  terminology  Is  justified  on  Intuitive  grounds 
as  well  as  on  the  grounds  that  It  can  be  proven  that  If  a  charge  density  Is 
moment  expanded,  then  the  corresponding  expansion  of  the  electro- static 
potential  Is  the  conventional  multi  pole  expansion.  It  will  be  clear  that  It 
Is  necessary  to  moment  expand  the  two- point  tensors  if  one- point  scale 
equations  are  to  hold.  Me  observe  that  for  any  function  f(x)  of  the  vector 
variable  x^  the  Dirac  identity  holds: 

f(x)  =  J six  -  x')  f  (x')  dx*  (67) 

6  Is  even  In  Its  vector  argument  and  by  Schwartz's  theory  it  can  be  Taylor 
expanded.  According  to  Schwartz  distribution  theory,  the  convergence  of  the 
following  expressions  Is  understood  as  convergence  In  measure.  We  can  thus 
write  the  two  expansions: 


6(x  -  x')  *  6(x)  -  x}  —  6< x)  +  j  xjxi  ■  6(x)  +  ... 

'  ax.j  c  1  j  ax^axj  — 


(68) 


and 
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=  6(x')  -  xi  -rr  6(xi)  +  -  XjXj  — T7~T  6( x 1 
—  1  ax^  -1  2  1  j  ax}axj  - 


6(  x  * )  + 


(69) 


Substitution  of  Eq.  (68)  Into  (67)  yields 


f(x)  *  a(x) 


/■ 


f(x')  dx' 


f(x')  dx* 


+  ... 


(70) 


which  we  call  the  moment  expansion  of  f.  Substitution  of  Eq.  (69)  Into 
Eq.  (67)  yields,  on  the  other  hand. 


f(x)  *  f(0)  -  xi 


/jL\  + 

[dxj**0  •** 


(71) 


which  Is  the  standard  Taylor  expansion  of  f(xh  The  moment  expansion,  Eq. 
(70),  is  suitable  for  the  highly  peaked  function  while  the  Taylor  expansion, 
Eq.  (71),  Is  suited  to  the  opposite  regime  of  slowly  varying  functions. 

A  remarkable  and  useful  duality  can  be  proven  for  the  Fourier 
transformation  f ( k) : 


.  C  ik*x  ,  „ 

?( k)  =  I  e  “  -  f ( x)  dx 
Insertion  of  Eq.  (70)  into  Eq.  (72)  yields 

f(k)  *Je--  6(x)Jf(x')  dx'  - 

*  I  f U*  )dx'  +1k^  lx}  f( jc*  )dx'  +  ... 


+  ... 


dx 


*  f(0)  +  k; 


where  we  have  used 


+  ... 


«0)  =  Jf(x')  dx' 
(^7-)  =  -  fxi  f(x)  dx' 

\aki/k*0  J 


which  can  be  deduced  from  Eq.  (72). 

Insertion  of  Eq.  (71)  into  Eq.  (72)  on  the  other  hand  yields 


f(k)  *  f e1-*-  f(0)  -  +  ... 

J  L  \dki/x=o 


dx  = 


(72) 


(73) 


(74) 


(75) 
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(2k)3  6(k)  f(0)  +  1(2k)3  (-^]  + 

“  "  a*i  \axi/x*0 

k)Jf(k')  dk*  -  ?  (k‘>  dk  +  ... 


where  we  have  used 


f(0)  =Jf(u')  dk7(2n)3 


(76) 


(77) 


Jk\  f(k')  dk'/(2n)3  (78) 

which  follow  easily  from  the  inverse  of  Eq.  (72). 

We  may  now  expand  the  tensors  R^j,  R^,  and  R  in  terms  of  their  moments. 
From  dimensional  considerations  the  moments  necessarily  define  lengths.  A 
simple  choice  was  suggested  by  Rotta  (Ref.  2)  for  the  velocity  correlation 
tensor.  This  choice  was  alluded  to  earlier  in  this  paper  and  is 


where,  of  course 

q2(xc)  *  A*1"  Rkk^-c*  t)  (31) 

r+0 

An  alternative  Identification,  convenient  for  certain  purposes  (i.e.,  the 
homogeneous  flows  of  Section  2),  is 
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(82) 


which  has  the  same  dimensionality  but  a  somewhat  different  tensor  character 
than  the  Rotta  version,  Eq.  (80).  We  adopt  In  analogy  with  Eq.  (78) 
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The  analogues  of  the  identification,  Eq.  (80),  for  the  velocity  scale 
may  be  chosen  to  be 


(1)  , 
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and 
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(86) 


The  leading  term  in  the  moment  expansion,  Eq.  (70),  can  then  be  taken  to  be 

(2) 

R1j(Xc,r)  =  Mi j  (x^  6(r)  (87) 

where  we  have  Introduced  the  one-dimensional  Dirac  function  6(r)  (r  is 
defined  as  the  magnitude  of  the  vector  r_,  r  *  r).  The  main  property  of 
fi(r),  for  the  present  purposes,  is  its  relation  to  &(£)•  This  is  given  by 


6  ir) 
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(88) 


We  also  Introduce  the  moment  expansion  for  the  correlation  Rj  and  R  as 


fol 1 ows 
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(89) 
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and 


(0) 

R(xc,  r)  a  M  (xc)  6(r)  (90) 

In  order  to  obtain  the  rate  equations  for  the  moments,  we  require  three 
steps  which  we  are  now  In  a  position  to  take: 

(a)  Develop  a  closure  model  for  the  R^j  equation  so  that  when  ^  -►  x, 
we  obtain  a  standard  equation  for  u}uj.  This  Is  a  useful  correspondence 
criterion  designed  to  allow  us  to  Incorporate  what  has  already  been  gained 
from  second-order  closure.  Of  course,  what  Is  said  for  R^j  is  taken  to 
apply  to  R^  and  to  R  as  well. 

(b)  Substitute  into  the  closed  equations  for  the  two-point  correlation 
tensors,  R^j,  R^  and  R,  their  moment  expansions  and  retain  the  leading  term 
to  define  a  first  approximation  to  the  behavior  of  these  quantities.  We 
observe  that  the  moment  expansion  appears  here  as  an  intermediate  step, 
prior  to  Integration  over  the  variable  £w1th  respect  to  which  we  wish  to 
average.  We  return  to  this  essential  point  below  (c). 

(c)  We  Integrate  over  all  jr  the  first-order  expanded  model  equations 
for  the  correlation  tensors.  This  step  of  integration  (or  averaging)  is 
taken  In  many  known  moment  procedures  without  the  intermediate  step  (b).  We 
return  below  to  this  point  to  clarify  our  position. 

At  this  juncture  we  believe  that  two  main  observations  are  worth 
making.  With  respect  to  (a),  we  observe  that  in  order  to  have  a  closed  set 
of  partial  differential  equations  which  are  supplemented  only  by  boundary 
and  Initial  conditions  we  must  obtain  closure  at  the  level  of  the  two-point 
correlation  tensors.  This  Is  true  provided  we  identify  the  spatial  scales 
Introduced  In  the  (one-point)  second-order  closure  equations  as  integral 
scales  associated  with  two- point  correlation  functions  when  we  utilize  the 
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correspondence  principle  expressed  in  (a).  We  are  cognizant  that  such  a 
procedure  is  a  strong  assumption.  Hence,  our  label:  "generalized 
second-order  closure." 

The  second  point  that  we  wish  to  make  refers  to  the  requirement  of  an 
intermediate  step,  the  moment  expansion  in  (b),  prior  to  the  development  of 
the  moment  equations.  This  step  is  the  least  orthodox  in  our  development, 
and  we  have  repeatedly  attempted  to  bypass  it.  We  believe  that  the  moment 
expansion  is  necessary  for  the  development  of  rate  equations  for  the  scales. 
Reconsider  briefly,  as  a  prototype,  the  standard  procedure,  due  to  Maxwell, 
utilized  in  obtaining  moment  equations  from  the  Boltzmann  equations  (for  the 
single-particle  velocity  distribution  function  in  a  gas)  on  the  way  to 
deriving  the  transport  coefficients  for  the  Navier-Stokes  equations.  We 
consider  here  velocity  space  moments  of  the  quantity 

6f  =  f(v,t)  +  v»yf  -  Jtf,f]  (91) 

at  —  — 

where  J  is  the  collision  integral,  bilinear  in  f.  That  is,  we  multiply 
Eq.  (91)  by  powers  of  \[  and  then  integrate  over  all  velocity  space.  We 
observe  that:  (i)  the  convection  term  is  linear  in  the  averaged  variable  v^ 
and  hence  raises  the  moment  index  by  precisely  one,  and  (ii)  the  first  five 
moments  of  J  vanish  as  a  consequence  of  conservation  of  mass,  momentum  and 
kinetic  energy.  An  orderly  hierarchy  of  moments  ensues,  with  the  moment 
equations  for  momentum  and  for  energy  containing  exactly  one  higher  moment 
each.  The  closure  of  these  equations  and  hence  the  calculation  of  these  two 
higher  moments  is  then  the  task  of  the  Chapman-Enskog  expansion.  Feature 
(i)  is  characteristic  of  all  the  moment  procedures  known  to  the  writers,  and 
it  is  essential  for  the  implementation  of  Maxwell's  approach.  The  absence 
of  this  feature  in  the  equations  for  the  correlation  tensors  of  turbulent 
fluids  is  what  forces  our  intermediate  step,  the  moment  expansion.  For  the 
correlation  tensors,  convection  occurs  independently  at  two  distinct  points. 
By  contrast,  the  integral  scales  evolve  at  one  (and  each)  point  of  space 
since  they  are  functions  of  the  single  variable  Xc*  We  submit  that  the 
reduction  of  the  two  points  x  and  £  to  the  single  x^.  cannot  be  accomplished 


without  our  moment  expansion,  Eq.  (70). 
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The  reader  can  rather  easily  convince  himself,  in  the  cases  where  the 
Intermediate  step  (b)  of  carrying  out  the  moment  expansion  is  not  necessary 
(or  in  the  case  described  above  of  Maxwell  moments  of  the  Boltzmann 
equation),  that  the  explicit  use  of  step  (b)  causes  no  deviations  at  all 
from  the  standard  procedures. 

When  we  carry  out  the  program  contained  in  (a),  (b)  and  (c)  above,  we 
obtain  the  sought-after  rate  equations  for  the  moments  (at  the  single  point 
x^.).  We  illustrate  the  procedure  with  the  simplest  term,  the  partial  time 
derivative  of  R^j.  Application  of  the  procedure  to  other  terms  in  the 
two-point  correlation  equation  is  tedious  but  not  difficult. 

(a)  Consider  the  relevant  term  in  the  modeled  equation.  The  term  is 


JL  r..  (x  r) 
ij  "at  u  - 


(92) 


For  this  term,  of  course,  no  modeling  was  necessary. 


(b)  Insert  in  the  term  (modeled  or  not)  the  moment  expansion  and 
retain  the  leading  term 
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where  we  have  used  Eq.  (87).  It  may  not  be  amiss  to  remark  that  a  full 
series  expansion  can  be  considered  even  though  for  simplicity  we  focus  here 
on  the  leading  term. 

(c)  Integrate  the  quantity  obtained  in  (b)  over  all  space  of  relative 
position  r  using  a  weight  appropriate  to  the  generalized  moment  considered. 
Using  the  Rotta  moment  defined  in  Eq.  (79),  we  consider 
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i  used  the  basic  property,  Eq.  (88),  of  6( r)  and  the  familiar 


djr  s(jr)  =  1 


(95) 


see  that  implementing  steps  (a),  (b)  and  (c)  yields  rate 
•  the  moments. 


general  formula  can  be  given  for  averaging  with  a  function  A 
triable  which  can  be  the  component  of  a  tensor 
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der  as  a  complete  example  the  following  model  equations  in 
mplicity,  diffusion  terms  have  been  omitted.  We  thus  consider 
neous  flows  at  high  turbulent  Reynolds  number.  For  the 
ocity  correlation  tensor  we  assume 

(aui  au<  \  n 

'* ,j]  + 
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For  the  two-point  flux  of  the  passive  additive  we  assume 


(97) 


=  -  Rik  r1-  -  —  Rk  +  v'  ?  tRi  -  M  -  A  ?  CR,-  -  Tj  (98) 
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and  for  the  two-point  autocorrelation  function  of  the  passive  additive 
(e.g.,  temperature  autocorrelation),  we  assume 


fii , .  [r  -aa.  t  j£  R 1 . 
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where  the  two-point  material  derivative  D/Ot  is  defined  by 


JL^  =  —  +  u  (x)  — +  u.,(y)  — — 
^Dt/two-  at  k  -  ax  kv^  ayk 
point  k 


(100) 


The  model  parameters  v,  v' ,  A,  b  and  s  will  be  discussed  in  Section  2.  In 
addition  to  Eqs.  (97),  (98)  and  (99),  we  consider  the  requirements  of 
continuity: 
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We  then  obtain,  from  the  two-point  models  given,  the  following  one-point 
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correlation  equations  by  taking  the  limit  r_ 
stress  <uiui>  satisfies 


The  kinematic  Reynolds 


The  flux  of  passive  additive  satisfies 


and  finally  the  temperature  variance  satisfies 


—  e1^  =  -  2  uTe1 
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We  have  used,  in  the  one-point  correlation  equations,  Eqs.  (104),  (105)  and 
(106),  the  one-point  material  derivative 


te) 


one- 

point 


=  £  *  uk  <*c> 


(107) 


omitting  the  distinction  between  Eqs.  (100)  and  (107)  since  no  confusion 
can  arise.  The  equations  (104)  through  (106)  are  those  of  our  standard 
model  as  given  by  Lewellen^. 

When  we  apply  to  the  two-point  equations,  (97),  (98),  and  (99),  the 
moment  expansion  and  average  procedure  specified  above,  we  obtain  the  moment 
equations.  A  tensor  scale  is  contained  in  the  velocity  moment  equation 
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A  vector  scale  Is  determined  by  the  moment  of  the  flux  which  satisfies 
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And  finally  the  temperature  scale  Is  determined  by 
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We  shall  see  In  Section  2  that,  accepting  the  Identification  of  the  scales 
given  In  Eq.  (80),  we  obtain  good  agreement  with  experiments  on  near- 
homogeneous  turbulent  shear  flows. 


2.  NEARLY  HOMOGENEOUS  FLOWS 


The  goal  of  this  section  is  to  demonstrate  that  the  theory  of  turbulent 
scales  proposed  in  the  previous  section  is  in  considerable  agreement  with 
the  experiments  carried  out  at  JHU/APL  on  nearly  homogeneous  flows. 

To  introduce  our  discussion  of  nearly  homogeneous  flows,  we  observe 
some  general  features  of  the  scale  equations  that  we  have  derived  and  then 
will  exhibit  exact  solutions  for  theoretically  homogeneous  flows.  Adopting 
our  definition  of  the  tensor  integral  scale  for  velocity  components  that  we 
discussed  in  the  previous  section,  namely. 
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the  moment  equation  for  M:.' 
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can  then  be  written  as  (we  drop  the  subscript 


(112) 


As  indicated  in  the  Introduction,  it  is  promising  to  identify  the  scale  A 
that  enters  the  modeled  terms  with 

A  =  J%  (113) 

2 

We  can  contract  Eq.  (112)  and  subtract  the  rate  of  q  to  obtain  a  local 
rate  equation  for  A  .  This  reads 
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It  is  of  interest  that  the  coefficient  of  3uj/3Xj  in  the  first  term  on  the 
right-hand  side  of  Eq.  (114)  is  not  proportional  to  uTuT  .  Such  as  assump¬ 
tion  has  been  made  repeatedly  in  the  universal  "scalar  scale"  models.  A 
simple  approximation,  that  we  call  the  a  approximation,  can  be  deduced  from 
the  assumption 
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Substituting  the  ansatz  Eq.  (115)  into  Eq.  (114),  the  "production  term" 
simplifies,  and  we  find 
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Furthermore,  we  can  obtain  an  equation  for  a  by  substituting  the  ansatz 
Eq.  (115)  into  (112),  subtracting  the  rate  of  change  of  A  given  by  Eq. 
(116)  and  contracting  the  resulting  second-rank  tensor  equation  by  the 
deviator  <(u!uj)>-  (l/3)6ij.q2  and  solving  for  3o/3t  .  The  resulting  a 
equation  is 
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The  absence  of  nondiffusive  equilibrium  for  A  as  given  in  Eq.  (116) 
is  by  no  means  accidental.  It  is,  in  fact,  essential  for  a  good  match  of 
the  model  parameters  to  the  nearly  homogeneous  flow  data  (both  for  simple 
grid  and  for  near-homogeneous  shears).  We  note  from  the  a  equation,  as 

O 

pointed  out  by  Donaldson  and  Sandri  ,  that  flows  with  different  turbulent 
structure  have  a  different  coefficient  (a/3)  -  1  in  front  of  the  A  scale 
"production"  term.  Thus  a  scalar  scale  equation  which  does  not  take  into 
account  the  turbulent  structure  cannot  be  valid.  Furthermore,  we  can  prove 
that 

4 

F  =  <u!ul><u!up  -  ^->0  (118) 


quite  generally,  with 

q2  -  <u'u'>  (119) 

Therefore,  the  only  non-negative  stationary  solution  of  Eq.  (117)  which  is 
nondiffusive  is  a  =  3  .  We  shall  see  below  that  all  the  solutions  for  homo¬ 
geneous  shear  asymptote  to  a  solution  that  can  be  exhibited  exactly  and  for 
which  a=3.  It  is  interesting  that  the  value  of  a  used  in  practice 
(in  the  universal  scalar  scale  modeling)  is  a  =  2.48  .  As  pointed  out 

O 

again  by  Donaldson  and  Sandri  ,  it  is  attractive  to  conjecture  that  diffusion 
always  lowers  the  value  of  a  from  its  asymptotic  value  for  homogeneous 
flows  and  that  the  different  behavior  of  two-dimensional  and  axisymmetric 
jets  can  be  attributed  to  the  substantial  differences  in  the  diffusion  terms 
in  the  two  geometries. 

We  now  present  two  exact  solutions  for  ideally  homogeneous  shears.  We 
start  with  some  notation. 
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The  centroid  vector  and  mean  velocity  vector  are  taken  to  have 
components 

(x,y,z)  ,  (U(y),0,0) 

with  3U/3y  =  U'  =  constant.  The  relevant  components  of  the  Reynolds  stress 
equations  are  obtained  from  Eq.  (46).  We  drop  primes  on  the  fluctuations  and 
give  a  form  useful  for  numerical  integration  in  which  u^  and  Ajj  are 
calculated  from 

~T  T  ~T  ~T 


U1  =  ^  '  u2  -  u3 


(120) 


Al!  =  3A  '  A22  '  A33 


(121) 


The  other  relevant  components  of  the  stress  and  energy  equations  are 


it  *?"  j (1  ■ 2b)  V" 


ltT=I(1'2b)  Wu3 


(122) 


(123) 


3t  U1U2  'u2  U'  "  A  U1U2 


(124) 


aa..  Viu-  .bd 

3t  q  u  A 


(125) 


For  the  tensor  scale  components,  we  obtain,  using  Eq.  (112), 


3t  “22 


\  A22  +  (1  -  2b)q 


(126) 


Aoo  =  ■  r  A-31  +  (1  ~  2b)q 


(127) 


(128) 


3t  A12  =  '  T  A12  "  U'A22 


3  ulu2  2 

ft  A  =  2^U,A.|ai2U'  +  Vq 


(129) 


where 


1  u,u,, 

j  =  -2  ^-U'  +  J  (1  -  2b  -  v1) 

q 


(130) 


2.1  Solution  of  the  Shearless  Equations 

Setting  U'  =0  ,  we  see  that  equations  for  q  and  A  decouple  from 
the  tensor  components.  Introducing  the  deviators 


1 


d- .  =  u .u .  -  4-  6 .  .q‘ 
ij  i  J  3 


(131) 


Dij  Aij  ~  3  6ijAkk 


(132) 


and  the  time 


-I 


(133) 


we  have  the  set 


3t  q  ~  "  T  q 


ft A 


(134) 


3 j  _  b  *  3  n  _ 

3t  ij  '  “  x  ij  ’  3t  ij  ~ 


1  -2b-v' 


(135) 


ft  T  =  b  +  v‘ 


(136) 


Integrating  Eq.  (136)  as 


x  =  (b  +  v')(t-  tj  + 


o'  q 


(137) 


we  see  that  q  ,  A  ,  d^  and  D_  are  suitable  powers  of  (q0/T0)x  ;  for 
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example. 


J 


1 


i» ■ 

1 


q  =  q. 


d. .  =  d. . 

ij  ij, 


A  =  A. 


-i-b/(b+v' ) 


(b  +  v* )  (t  -  tQ)  +  1 
o 


-i  -l/(b+v* ) 


(b  +  v*)  ^2.  (t  -  tQ)  +  1 
o 


-T-v'/(b+v' ) 


(b  +  V1)  /  (t  -  t0)  +  1 
0 


(138) 


(139) 


(140) 


A  good  fit  to  experimental  data  on  the  decay  of  turbulent  energy  and 
the  growth  of  eddy  size  for  grid  turbulence  is  obtained  if  one  chooses 


b 


(141) 


v'  =  0.075 


We  then  see  that  for  large  times 


(142) 


(143) 


which  shows  that  the  deviator  decays  with  a  power  about  four  times  larger 
than  the  energy. 

From  the  solutions  given  above,  we  can  verify  that  statistics  are  pre¬ 
served  by  the  model  equations  if  the  model  parameters  satisfy  certain  bounds. 
We  first  show  that  the  two  tensors  uTiiT  and  A.,  are  positive  definite 

*  J  *  J 

from  their  definitions.  Consider  an  arbitrary  (constant)  A.,  then, 

AiVTA.Mu.A)2>0  (144) 

the  equality  sign  holding  for  A  =  0  only.  Thus,  uTuT  is  a  positive 
definite  tensor. 
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From  the  definition  of  the  scale  tensor  given  by  Eq.  (Ill),  using 


Fourier  transform  on  R..  , 

■  J 


.2  /*r..  /v. 

1- A  =  I  _KL  d;  =  I 

3  1J  J  4*7  J  8lTk 


3idt 


(145) 


where  the  power  spectrum  tensor  4>..  is  positive  definite  by  Khiutchine's 

'  J 

theorem.  Thus, 


\  Vi  jAj =  m  Ai<i>i  jAj  - 0 


(146) 


Thus,  is  positive  definite  because  q  is  positive  as  a  consequence 

of  Eq.  (144). 

Using  the  solution  (141)  and  an  analogous  solution  for  A..  ,  we  find 

■  J 

v  -l/{b+v<) 


Vj(t>  * 


r/o  \  -2b/(b+v* )  /  \  -l/(b+v' ) 

-ft-) 


(147) 


Aij(t)  -  Aij(0)  (  fQ' 


-(l-2b-v’ )/(b+v‘ ) 


■/.  \  v'/(b+v' )  /  \ -(l-2b-v‘  )/(b+v' )" 

•"'•'.i  ft’)  -ft’  j 


(148) 


We  now  multiply  Eq.  (147)  by  when  A^  is  an  arbitrary  vector,  and 

find 


Substituting  these  forms  into  the  differential  equations, 
we  find  that  the  exponentials  cancel  and  that  an  algebraic  set  of  equations 
for  the  amplitudes  is  obtained.  It  is  possible,  with  some  algebra,  to 
solve  the  amplitude  equations  explicitly  in  terms  of  the  parameters  b  and 
v1  .  The  energy  components  are 


7 

~T~ 

q 


1  +  6v*  +  4b 
3(1  +  2v  * ) 


=  0.5652 


(154) 


_2 

2 


3 

2 


1  -  2b 


1 


1  +  2v* 


=  0.2174 


(155) 


The  scale  components  are 


=  l±ivL±jb  .  1>696 

A  1  +  2v‘ 


(156) 


A22  .  A33  1  -  2b 

A  "  A  1  +  2v‘ 


0.652 


(157) 


We  see  that 


~T 

A11  _  U1  _  1  +  6v'  +  4b 
A 1  -  2b 


22 


T 

u2 


=  2.60 


(158) 


The  off-diagonal  components  are 

1 


Br  = 


!  ujU2 i 


1  +  2v 1 


(1  -  2b) (b  +  y')  _ 


0.194 


(159) 


1r = '  ttW  ^ 3(1  ■  2b)(b  +  v')  =  -o. 


194 


(160) 


The  Corrsin  parameter  is 
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/  1  +  4b  +  6v 


K 


2  ~2 
U1  u2 


The  ratio  of  the  two  times  is 

1  =  _5_=  _J _ 

a  U'A  1  +  2v‘ 


/I  -  2b 
3 ( b  +  v’ 


0.972 


(162) 


and  the  growth  rate,  a  ,  is 


,  _  w.  I  _  1  /I  -  2b 

a  "  v  a  ~  1  +  2 v 1  /3fb  +  v' 


=  0.073 


(163) 


We  notice  two  additional  interesting  parameters: 


I  ulu2  *  u1 

— =  ot  *  Br  =  b  +  v'  *  0.20 


(164) 


hi  X  =  3(  _  1  5 

A  6ij  \  q2  3  6ij 


(165) 


In  terms  of  the  o  parameter,  Eq.  (165)  corresponds  exactly  to  the  value 


o  =  3 


(166) 


As  remarked  above,  o  =  2.48  has  given  satisfactory  results  on  a  number  of 

flows.  Also  note  that  for  this  coupled,  asymptotic,  convective  solution, 

2 

the  ratios  <u.u.>/q  ,  q/All '  are  modified  from  those  of  the  equilibrii 

1  J  g 

solution  with  no  diffusion  (the  superequilibrium  introduced  by  Donaldson). 

Por  these  solutions,  we  have  (using  Eq.  (46)), 


JT  3 


(167) 


< vv>  _  <ww>  _  1  -  2b  _ 


< 


[ij 


q  =  1  -  2b  =  2 

A(9U/3y)  3b 


(17G) 


now  appropriate  that  we  consider  the  agreement  between  these 
results  just  obtained  and  experimental  measurements.  The  meas- 
will  make  use  of  are  those  of  Harris,  Graham  and  Corrsin  (Ref.10). 

re  is  Figure  3  from  Reference  10.  It  shows  the  growth  in  <^u.2^  , 

/  2  \  N  1  ' 
<^u2  /  in  a  constant  shear  flow.  For  these  measurements,  Uc 

/sec,  U'  =  3U/3y  =  48  sec"^  and  h  =  30.48  cm.  Also  shown  in 

is  the  behavior  of  the  Corrsin  number 

-  <ulu2^/ <u12Xu22> 

pect,  after  an  initial  transient,  that  the  flow  adjusts  to  an 
solution  with  certain  parameters  constant,  we  plot  faired  results 
3  in  semi -logarithmic  form  in  Figure  4.  It  is  apparent  from 
at  an  exponential  solution  has  been  reached  at  an  x/h  of  approx- 
We  note  that  the  growth  rates  in  terms  of  x/h  may  be  written 

<u2>  ,  <u|>  ,<u32>«e  c  (171) 


re  we  have  indicated  the  ac's  that  best  fit  each  curve.  An 
these  values  gives  a  =  0.087  or,  let  us  say,  a  =  0.09  .  To 
se  results  with  the  theoretical  results  just  given,  we  must  multi- 
Uc/hU’  =  0.848  which  gives  a  =  0.076.  This  is  not  in  bad 
ith  the  theoretical  value  of  0.0729. 

'am  has  been  written  to  solve  the  full  set  of  coupled  <u.u.>  and 
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>ns,  namely,  Eqs.  (46)  and  (112),  and  the  result  of  a  computa- 
versus  both  dimensional  time  and  nondimensional  time 


U1  At  = 


(172) 


Figure  5.  It  is  seen  that, although  the  growth  rate  of  q  is 
;ed,the  actual  values  of  q  are  some  20%  in  error  and  remain  so 
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in  the  asymptotic  region  of  the  solution.  This  is  not  good,  but  we  will 
discuss  a  possible  cause  of  this  error  after  we  have  exhibited  the  agreement 
between  experimental  results  and  all  the  parameters  we  have  derived  for  the 
exponential  behavior  of  a  homogeneous  shear  flow.  These  comparisons  are 
shown  in  the  second,  third  and  fourth  columns  of  Table  1. 

It  is  not  difficult  to  show  that  other  simpler,  second-order  models  of 
turbulent  flow  have  exponential  asymptotic  solutions  in  the  case  of  homo¬ 
geneous  shear  flow.  In  Table  1,  we  also  show  the  theoretical  results  for 
the  following  models: 

(a)  A  full  closure  in  the  case  of  the  <^u^>  tensor  and  the  single  scale 
equation 

^  =  0.35  <uv>  U'  +  v 1  q  (173) 

(b)  A  q  -  A  model  constructed  from 

2 

da  =  .  luv >..  u-  .  bg_  fl74\ 

dt  q  A  *1/4' 

^  =  0.35  4%  <  uv>  U'  +  v'q  (175) 

q 

together  with  the  assumption  that 

<uv>  =  -0.35qAU'  (176) 

An  examination  of  Table  1  is  instructive.  First  of  all,  the  tensor 
scale  model,  with  only  two  adjustable  constants  which  were  set  from  an  exper¬ 
iment  on  grid  turbulence,  seems  to  do  the  best  job  overall.  It  gives  by  far 
the  best  prediction  of  asymptotic  growth  rate.  The  single  scale  model  as 
normally  used  (which  has  the  added  parameter  c  =  0.35)  is  pretty  good  on  all 
quantities  except  q/AU'  and  the  growth  rate,  where  it  is  very  poor  indeed. 
The  approximate  q  -  A  which  has  still  another  adjustable  parameter  does 
least  well  for  the  parameters  chosen,  which  are  typical  of  these  boundary- 
layer-  like  flows. 


A  most  interesting  result  is  found  if  we  make  the  assumption  that,  as 
indicated  by  the  theoretical  developments  we  have  presented,  the  "constant" 
in  the  simple  scale  equation  should  be  zero  and  not  equal  to  0.35  for  homo¬ 
geneous  shear  flows.  If  we  do  assume  c  *  0  in  this  formulation  and  find 
the  asymptotic  values  of  the  basic  parameters  fcr  homogeneous  shear  flow,  we 
obtain  almost  the  same  values  of  the  parameters  that  were  obtained  using  the 
tensor  scale  equation.  This  is  not  surprising  since  the  two  formulations  are 
now  very  similar.  However,  it  must  be  remembered  that  it  was  the  tensor 
scale  equation  which,  when  solved  for  homogeneous  shear  flow,  showed  that  the 
production  term  should  disappear  from  the  scale  equation. 

If  one  puts  c  =  0  in  the  q  -  A  model  we  have  concocted,  the  growth 
rate  and  the  Bradshaw  number  become  too  large  (0.963  and  0.257,  respectively) 
while  the  parameter  q/AlT  drops  to  1.284  (which  is  not  a  great  improvement). 


3.  CONCLUSIONS 


In  the  previous  sections,  we  have  reviewed  some  of  the  characteristics 
of  second-order  modeling  as  it  is  currently  used.  One  of  the  primary 
criticisms  of  these  methods  has  been  that  they  take  no  account  of  the 
structure  that  can  be  found  in  turbulent  fields  by  modern  instrumentation. 

We  have  given  here  an  outline  of  how,  by  the  use  of  a  simple  definition  of 
tensor  scale,  second-order-closure  models  might  be  extended  to  take  account 
of  information  on  structure  that  can  be  gleaned  from  the  two-point  correla¬ 
tion  equations.  The  tensor  scale  used  is  certainly  not  ideal  for  this  pur¬ 
pose,  but  is  was  used  not  only  because  we  are  familiar  with  it  but  also 
because  it  illustrates  many  features  that  will  be  exhibited  by  any  other 
definition  of  tensor  scale. 

We  believe  we  have  demonstrated  two  things  in  the  results  presented. 
First,  we  believe  we  have  shown  that  there  really  cannot  be  such  a  thing  as 
a  universal  scalar  scale  equation.  Hence  we  believe  that  any  steps  taken 
to  improve  second-order-closure  methods  in  the  future  must  include,  among 
other  things,  the  derivation (from  appropriate  models  of  the  two-point  corre¬ 
lation  tensor  equations)multiple  scales  which  will  give  a  model  that  is 
compatible  with  the  structure  of  the  turbulent  eddies  that  exist  in  a  given 
mean  flow.  The  method  we  have  used  here  defines  a  tensor  scale  and  uses  a 
moment  expansion  to  look  at  some  general  features  of  the  structure  problem 
that  can  be  derived  from  a  particular  definition  of  tensor  scale.  The  method 
is  a  good  approximation  for  homogeneous  flows.  It  is  less  justifiable  for 
nonhomogeneous  flows.  However,  we  believe  that,  at  the  present  time,  it 
bears  a  relation  to  a  more  complete  formulation,  much  like  eddy  viscosity 
methods  bear  to  more  complete  formulations  for  calculating  the  Reynolds 
stress  correlation  uTuT  . 
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Second,  we  believe  that  we  have  shown  that  the  homogeneous  shear 
experiments  are  very  powerful  tools  for  the  modeler.  We  believe  that  they 
do  indeed  have  asymptotic  solutions  that  are  exponential  and  that  when  the 
asymptotic  state  is  reached  certain  nondimensional  parameters  become  constant. 
Since  the  asymptotic  value  of  these  parameters  can  be  computed  from  a  given 


model,  the  experimental  results  are  an  extremely  useful  tool  for  the  develop¬ 
ment  of  valid  models.  One  reason  the  measured  values  of  these  parameters  are 
so  useful  is  that  they  are  independent  of  initial  conditions  and,  in  the  past 
arguments  over  initial  conditions  have  been  used  to  cover  a  multitude  of 
modeling  sins. 


4.  DISCUSSION 


There  is  more  than  one  way  to  look  at  any  physical  problem.  In  the 
field  of  turbulence  research,  it  has  become  fashionable  to  define  almost  any 
unsteady  flow  in  the  wake  of  an  unaccelerated  body  as  turbulence.  A  lot  of 
these  unsteady  flows  contain  large-scale,  coherent  structures  and,  thus, 
there  has  arisen  a  "large-scale  eddy  cult."  The  basic  thesis  of  this  group 
is  that  turbulence  should  or  (for  the  extremists)  must  be  attacked  by  some 
technique  which  identifies  the  turbulence  with  the  interaction  and  decay  of 
such  large-scale  structures.  A  corollary  of  this  position  is  that  closure 
methods  do  not  or  (for  the  extremists  )  cannot  address  themselves  specifi¬ 
cally  to  large-scale  eddies  and,  therefore,  are  not  really  anything  but  dull, 
unphysical,  and  temporary  methods  for  dealing  with  turbulent  flows. 

The  authors  of  this  paper  do  not  believe  that  these  people  really 
understand  the  nature  of  closure  calculations  at  the  present  time.  Not  only 
have  closure  methods  demonstrated  the  existence  of  large-scale  eddies  in  two 
cases  (the  roll  eddies  of  the  marine  planetary  boundary  layer11  and  unsteady 
large  eddies  in  the  flow  behind  a  rearward  facing  stepAt),  but  these  eddies 
have  been  resolved  in  all  their  gory  detail.  This  can  be  accomplished  when 
the  closure  equations  are  used  in  their  elliptic,  time-dependent  form  and 
the  grid  spacing  used  is  fine  enough.  Why  is  this  so?  It  is  because  the 
Euler  equations  (which  govern  the  formulation  and  a  great  deal  of  the 
behavior  of  large  eddies)  are  contained  in  the  time-dependent,  elliptic 
equations.  Thus  we  submit  that  closure  techniques  not  only  can  describe 
large-scale  eddies  but  must  do  so  if  the  time-dependent,  elliptic  forms  of 
the  equations  are  used  and  the  grid  spacing  used  is  small  enough  to  resolve 
these  eddies. 

One  further  point.  The  Karman  vortex  wake  is  an  unsteady  flow  assoc¬ 
iated  with  an  unaccelerated  body.  We  do  not  prefer  to  think  of  it  as  a  form 
of  structured  turbulence  (nor  did  those  who  first  studied  the  phenomenon) 
although,  if  one  were  a  member  of  "the  cult,"  this  point  of  view  might  be 
taken.  The  reason  that  we  choose  not  to  consider  it  as  structured  turbulence 
is  that  there  is  ample  evidence  of  its  existence  in  laminar  flow.  From  the 


closure  point  of  view,  the  Karmen  wake  Is  looked  at  as  an  unsteady  flow 
peculiar  to  the  body  that  produces  It.  This  unsteady  flow  Interacts  with 
Itself  to  decay  either  through  laminar  exchange  and  dissipation  or  through 
turbulent  exchange  and  dissipation,  depending  on  the  Reynolds  number. 

The  research  reported  here  Is  a  description  of  our  first  attempts  to 
make  a  closure  theory  of  turbulence  that  Is  compatible  with  the  large-scale 
structures  that  we  must  Inevitably  find  when  we  run  our  closure  codes  In  an 
elliptic  manner.  The  work  reported  Is  one  completed  step  in  this  direction. 
As  noted  in  the  text,  we  do  not  consider  this  work  complete.  We  sincerely 
regret  that  we  had  to  terminate  this  work  before  our  attempt  to  construct  a 
more  general  closure  formulation  could  be  completed.  While  we  consider  many 
of  the  detailed  large-eddy  studies  to  be  outstanding,  we  also  believe  that 
elliptic,  nonsteady,  closure  techniques  will  be  the  backbone  of  turbulence 
computations  that  will  be  of  use  to  the  military  for  the  next  twenty  years, 
and  that  closure  techniques  that  are  compatible  with  this  approach  should 
and  will  be  pursued. 
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